ON THE LYAPUNOV EXPONENT OF A MULTIDIMENSIONAL 

STOCHASTIC FLOW 

MICHELE BALDINI f 

Abstract. Let Xt be a reversible and positive recurrent diffusion in R d described 
by 

t 

X t = x + ab(t)+ J m(X s )ds, 
'o 

where the diffusion coefficient a is a positive-definite matrix and the drift m is a 
smooth function. Let Xt(A) denote the image of a compact set A C R d under the 
stochastic flow generated by Xt- If the divergence of the drift is strictly negative, 
there exists a set of functions u such that 

lim / u(x)dx = a.s. 

t — foo J 

A characterization of the functions u is provided, as well as lower and upper 
bounds for the exponential rate of convergence. 

Keywords: diffusion, stochastic now, recurrence, super harmonic function, ellip- 
tic operator. 



1. Introduction 

In jHj Has'minskii proved that any solution of a one-dimensional stochastic differ- 
ential equation describing a positive recurrent diffusion Xt is stable with probability 
one in the metric given by the natural scales. More precisely, let Xt be a diffusion 
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with smooth coefficients 



t t 



(1.1) 



X t = x + / a(X s )db{s) + / m{X s )ds 







Then, for every pair of initial conditions y,z EM. 



(1.2) 



P lim " =0=1 



where s is the scale function. A different proof of this stability is offered in pQ, 
together with the result that the pathwise rate of convergence is exponentially fast, 
i.e. 



Here dll is the corresponding invariant measure. 

This paper proceeds in the same direction and tackles the multidimensional prob- 
lem. As expected, technical difficulties make the task quite challenging and more 
restrictive conditions on the coefficients are introduced. The main result is: Let Xt 
be the stochastic flow generated by a reversible and positive recurrent diffusion in 
M. d . If the diffusion coefficient is a positive-definite matrix and the drift is a C 1 func- 
tion subject to certain boundness conditions on the derivatives, then there exists a 
non-empty set of smooth functions u such that, for every compact set A C M. d , the 
u- volume of the flow Xt{A) tends almost surely to zero as t f oo. In addition, the 
upper and lower bounds for the corresponding Lyapunov exponent are provided. 

The paper is organized as follows: section 2 introduces the notation and the 
necessary background. Section 3 paves the road with some introductory results. 
The following two sections are devoted to the main theorems: section 4 studies the 
asymptotic behavior of the flow and section 5 discusses the Lyapunov exponent. An 
example is illustrated in section 6 and section 7 concludes. 



(1.3) 




on the lyapunov exponent 3 
2. Notation and preliminaries 

Throughout this article it is assumed that: d > 2 is the dimension of the phase 
space M. d , B is the cr-algebra of Borel sets of R d with the implicit assumption that 
every subset in R d belongs to B. Given D,D' C R d , D'ccD denotes a bounded 
subset such that D' C D. A is always a compact subset and Leb(^4) is the Lebesgue 
measure of A. 

Let D C M. d . A function / in D is called uniformly Holder continuous with 
exponent a G (0, 1] if 

x,y£D \ x y\ 

If the same inequality holds for every D' CC D, then / is called locally Holder 
continuous. A function / is called uniformly (locally) Lipschitz continuous if it is 
uniformly (locally) Holder continuous with exponent 1. Define the Holder space 
C k,a (D)(C k ' a (D)) as the space of functions whose A;th order partial derivatives are 
uniformly (locally) Holder continuous with exponent a G (0, 1] in D. 

Let (f2, J-, P) be the canonical Wiener space of the d-dimensional Brownian mo- 
tion with paths b(t), t > and let 6t be the measure-preserving shift operator on 
the space of Brownian paths, i.e. 9tb(s) := b(s + t). Let a be a d x d positive-definite 
matrix and a := aa l . Let m : W 1 — > R d be a drift coefficient such that: 

(i) m G C 1 ' 1 ^). 

(ii) There exists a potential tp : M. d —> (0, oo) that satisfies 
(2.2) m = -aVlnY>, 

< 2 - 3 > /£-'■ 

(iii) There exists a constant A c > such that 

(2.4) (V • m)(x) < -A c for all x G M d . 
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Roughly speaking, condition (hi) says that the drift resembles a restoring force a la 
Ornstein-Uhlenbeck, i.e. m(x) = —x. The reader will also notice that the potential 
is defined modulo a positive multiplicative constant, so to assure (|2.3|) it suffices to 
verify that ifj~ 2 is integrable. 

Let Xt(x,b) be the flow determined by the ltd stochastic differential equation 
(SDE) 

t 

(2.5) X t = x + ab(t) + J m(X s )ds. 

o 

Here "flow" means that Xt(x, •) is a function of a single Brownian motion, the 
same for every x, so that for every t > 0, Xt(x,-) is a 1-1 map from R rf into 
M. d . More precisely, after a null set in Q has been weeded out, the family of maps 
{Xt(x,b) : t > 0} defines a stochastic flow of global diffeomorphisms, i.e. for every 
Brownian path b 

(i) M. d 3 x i— ► Xt{x,b) is a diffeomorphism in R rf for all i > and the inverse 
map is continuous in t and smooth in x. 

(ii) Xt+ T (x, b) = Xt(-,0 T b) o X T (x, b) for all non-negative t, r and x £ 

For an exceptional presentation of this topic and for the proof that uniformly Lip- 
schitz continuous coefficients generate a flow of global diffeomorphisms, the reader 
should consult [1], section 4.7. 

The infinitesimal generator Q of the diffusion introduced above can now be worked 
out. For any f £ C 2 (R d ), 

(2.6) Gf = ^ 2 V ■ (V" 2 aV/) = 1 V • aVf + m ■ V/, 
whereas the formal adjoint operator is 

(2.7) £/ = iv • [V^ 2 aV(V 2 /)] = iv ■ aV/ - m ■ V/ - (V ■ m)/. 

Thanks to (|2.6|) - (|2.7() . the meaning of Q2.2JI and ()2.3() becomes clear: the first 
guarantees that Q is self-adjoint in L? ( ^tct ) (this is the reversibility condition in 
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[2] ) , the latter says that the invariant density measure of X t is indeed an invari- 
ant probability density. For general questions on the existence and uniqueness of 
solutions of SDE, we refer the reader to jSJ. 

An allied flow is now presented: Let X\[x, b) be the flow generated by the SDE 

t 

(2.8) xt-, + .m-j»w}» 



with reversed drift — m in place of m. Since the regularity conditions on the coeffi- 
cients are the same, the family of maps {x\(x, b) : t > 0} defines a stochastic flow of 
global diffeomorphisms of However, it will be shown that the ergodic property 
of the corresponding diffusions x\ is diametrically opposite to the one displayed by 
X t . 

The diffusion x\ is not the main subject of this paper, however it is intimately 
related to Xt (for more of this see pp). In the next section the ergodicities of both 
processes are studied, in conjunction with the harmonic property of Xf's infini- 
tesimal generator: it turns out that this operator enjoys a well-defined criticality 
condition which will be the key for the proof of the main theorems in the conclusive 
sections. 

For convenience of the reader, the generator and the adjoint operator are spelled 
out in full form: 

(2.9) Qtf = ^- 2 V ■ (V 2 aV/) = iv • dVf - m • V/, 

(2.10) £«/ = l -V ■ [V 2 aV(^- 2 /)] = \y ■ aVf + m ■ V/ + (V • m)f. 

3. Ergodic properties of X t and x\ 

A criterion of R. Pinsky can be applied to study recurrence and transience in 
terms of integral tests. 

Theorem 3.1 (0, pp. 258). Let L = \V ■ aV + aVQ • V onR d , d>2 satisfying 
the following conditions on every D CC M. d : aij{x) G C 1,a (D), Q(x) £ C 2,a (D) and 
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Yli j=i a ij{x)viVj > 0, for all v G R d — {0} and all x G D. Define 

(3.1) El (x) = K^) e 2QW, 

(3.2) E 2 (x) = _jf e 2 ^>, 

(a l (x)x,x) 

(3.3) E^r) = J E 1 (r6)d6, 

S d-i 

where d6 is Lebesgue measure on S"* -1 . 

(i) If 

oo 

(3.4) I r^^J—dr = oo, 

/ Mr) 

then the diffusion corresponding to L is recurrent. 

(ii) If 

oo 

< 3 - 5) / r "%p) dr<00 ' 

1 

on a subset of S^ 1 with positive Lebesgue measure, then the diffusion cor- 
responding to L is transient. 

Corollary 3.2. X t is positive recurrent and x\ is transient. 

Proof. Let \ m and Am be respectively the smallest and largest eigenvalue of a. Then 
< A m < Am < oo- A straightforward computation gives A m |x| 2 < (ax,x) < Am|^| 2 

I 1 2 I 2 1 

and 'jjj < (a~ 1 x,x) < Now, by virtue of the integrability of ip~ 2 , it follows 

that 

(3 - 6) I f sd -J^(re)de dr = 00 and 

oo 

(3.7) / r o/ — dr < oo for almost all 6 G S* -1 . 

i 
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A direct application of theorem 13. 1 1 establishes recurrence and transience for X% and 
Xf, respectively. What remains is to verify the positiveness of Xt. This follows from 
the observation that ip~ 2 is the invariant probability density of Xt (see discussion 
in chapter 4.9). □ 

Some spectral properties of the operator £^ are explored next. By the transience 
property, both Q$ and £^ possess a Green function ([S], Th. 4.3.3). This implies 
that for each non-negative / G C^^W 1 ) which is non identically zero, there exist 
positive solutions u G C 2,1 (IR d ) of £$u = —f. The smallest such solution is given 
by uq(x) = L d G(x,y)f(y)dy where G is the Green function of £" ([5], Th. 4.3.8). 
Define the cone of all superharmonic functions as follows: 

(3.8) 5 = {0<ti6 C 2,1 (]R d ) :0>Ae C%(R d ), £ tt u non identically zero}. 

Denote (£ the principal eigenvalue of i.e. the supremum of the spectrum of the 
corresponding self-adjoint extension. Using Theorem 4.4.5 in |Hj and the transience 
property it follows that 

ott 7 , 

(3.9) <£ = inf sup (x) < 0. 

weC 2 ' 1 (R d ) ^eRd u 

Thanks to the negative divergence of the drift, one can prove more. 

Theorem 3.3. The principal eigenvalue <£ is strictly negative. 

Proof. As stated above, for each non-negative / G CQ(M d ) which is not identically 
zero there always exists a positive u G C 2,1 (M rf ) such that 2$u = — / < 0. Now, with 
g G C 2 (M d ), the following identity holds: 

(3.10) Spg 2 = 2g &g - g 2 V -m + Vg ■ aVg. 



Since the square root does not spoil the regularity of a positive function, apply Q3.1U[) 
to g = y/u obtaining 2y/u&y/u = £"ii + uV • m — V-y/u • aSJyfu. Using Q2.4JI and the 
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positive definiteness of a, we can write 

. A 



(3.11) 

The proof is completed observing that 
(3.12) 



< 

u ~ 2 



<S < sup — /=-(«) < -y 



□ 



4. Main results 

Denote X\ the i-th component of Xt and rrij the j-th derivative of the i-th com- 
ponent of m. Denote |-D-X$| := \DX t (x, b)\ the Jacobian of the map x i— > X t (x, b). 



Lemma 4.1. Let < u G C 2 

u(X t )\DX t \ = ti(x)exp 



Then 



' Vu(X T ) 1 } (Vu-aVu){X T ) A }{S$u){X T ) 
ado — — I „ , , r , dr - 



u(X T ) 



2 J u 2 (X T 
o 



dr 



Proof. From standard results [4] , the gradient flow satisfies q-* 
Simple differentiation gives 



dXb. 



8X} 



ax? 



k=l 



dx H dx %d 



dt 



(4.1) 



V-m(X t )|DX t |dt, 



dx %l dx ld 



dt 



where V is the set of all permutations (1, . . . , d) — > . . . , i^) and p is the sign of 
each permutation. Now use Ito's lemma on u(Xt)\DXt\ to show 



(4.2) d u(X t )\DX t \ 
= \DX t \ 



Vu(X t ) ■ adb + ( - V • aSJu + m ■ Vu + (V • m)u ) {X t )dt 
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The proof is completed observing that the operator acting on u on the RHS of (|4.2jl 
is the formal adjoint of QK □ 



Next we turn to the analysis of the volume of Xt(A), i.e. the image of a compact 
set J 4cl d under the flow X t (-,b). 

Theorem 4.2. Let < u G C 2 (R d ) such that £fiu < 0. Then, 
(4.3) Ju:= Ju( X ) dx 



X t (A) X t (A) 



is a positive supermartingale. 



Proof. Denote Mt := u(Xt)\DXt\. For every < s < t and n > 1 define 

d 



(4.4) 



S n = min inf < r > s 

i=l,...,d 



/c=i 



n(X r ) 



> n 



and integrate (|4.2I) on [s, t A S* n ], 

iA5„ 

(4.5) M tASn =M s + [ M T Vu{ * r) .adb + 



tAS n 



u(X T ) 



Since the stochastic integral has expectation zero and &u < 0, then E[M tASn |j;] < 
M s . Because S n | oo, an application of Fatou's lemma gives E[Mt| < 7 r s ] < M s . 

For almost every Brownian path, M t is a positive and continuous map in the 
variable x, hence its integral is well defined. By the formula for the change of 
variables we have $ Xt (A) u = Ia^ anc ^ then using Fubini's theorem 



(4.6) 



E 



E 



Xt(A) 



< I M s = J u. 

A X S {A) 



□ 



Doob's theorem guarantees that limj^oo Jx t {A) u ex i s ^ s almost surely. It is shown 
next that it is always possible to find many n's such that the limit is zero. Condition 
(|'2.4[) . as expected, will play a crucial role in the proof. 
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{0 < 


u G C 2 (R d ) : 


u = 


p 


lim / u = 


= 1 




a— >oo y 











ueW, 
(4.7) 



Proof. By theorem 13.31 £^u < — 4fu. Set s = in (|4.5|) and take the expectation to 
obtain 

tAS„ tASn 

(4.8) E(M tA5 J = M + E J M T (£ ^^ r) dr < M - ^ E J M T dr. 

o o 

Sending n | oo and using Fatou's lemma we have 

t 

(4.9) E(M t ) < M - yE /" M T dr. 

o 

Next integrate over A, 

t 

(4.10) E J Af t < y M - yE /" (y M T )dr, 

A A OA 

and observe that the second term on the RHS is monotone decreasing and bounded 
from below, hence E f^°(J A M T )dr < oo. Therefore L (£ 4 M r )dr < oo almost 
surely. Since J Xt (A) u = Ia^ ^ s m tegrable and convergent as t j oo, the limit is 
necessarily zero. □ 

5. Lyapunov exponent 

In this section the rate of convergence of Jx t (A) u * s ana ly ze d to derive lower 
(upper) bound for the lim inf (lim sup) of 

(5.1) j In y u. 

Xt{A) 

The lower bound is easy to obtain. 
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Theorem 5.1. Let < u € C 2 (IR d ) and M t := u{X t )\DX t \. If 

Vn • aVu 



(5.2) 



< oo and 



|£«n| 
— 72" < °°' 



n 



then 



(5.3) 



lim - In Mt 

t— >oo f 



1 /" Vn • aVn 



u 2 ip 2 



+ 



a.s. 



and 



a.s. 



5.4 limmf-ln / n > Leb(A) -- / rp ^ + / — 

t^co t J 2 J u A ip z J uip 

X t (A) L Rd Rd 

Proof. To begin with, recall from lemma ETT1 that 
(5-5) 

{ fVu(X T ) „ 1 ? (V«-aVu)(X T ) , ?(£*n)(X T ), . 

\0 

It is easy to show that the process 



(5.6) 



B 



Vu(X T ) 
u(X T ) 



■adb 



is a standard 1-dimensional Brownian motion run with a random clock t(t) = 
Jo ^ U ^x) Xt ^ t ' virtue of the ergodic theorem for recurrent diffusions [7j (see 
also UJ, Th. 4.9.5), lim^oo ^ = J Rd ^ V " 2 °J"' > almost surely, leading to 



(5.7) lim I In (m^A = -\ I + [ 



uip" 2 



a.s. 



To conclude the first half of the proof, it must be checked that almost surely 
limsup^^ = hminfj^oo ^j^- = 0. Since 



(5. 



B(t) B(t) /21nlnt t 



t \/2tlnlnt V t t 



and almost surely lim^oo ln nt = 0, an application of the Hincin's law of the iterated 
logarithm proves it. 
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To prove (|5.4j) first apply the reversed Jensen's inequality for concave functions 
to obtain In j.Mt > fAnMt. Then Fatou's lemma permits to let t f oo under the 
integral sign, producing 

Iff 1 

(5.9) liminf — In u> liminf — In Mt a.s. 

t-^oo t J J t^oo t 

X t (A) A 

Formula (|5.3j) supplies the rest of the proof. □ 
The following result identifies the upper bound. Condition (|2,4|) is critical. 

Theorem 5.2. Let < u G C 2 (R d ) such that £$u < -^fu, then 

If A 

(5.10) limsup — In / u < — — a.s. 

t-*oc t J 2 

Proof. First recall that theorem 13.31 guarantees the existence of infinitely many 
smooth u > such that &u < -4fu. Next define 3* := e Xct/2 M t and let S n 
be as in Section 0J Ito's Lemma gives 

tAS„ 

(5.H) 3tASn<3,+ / 3r^^-^d6, 

J u{X T ) 

s 

permitting to conclude, after an application of Fatou's lemma, that 3t and j A "bt 
are both positive supermartingales. It follows that there must exist a finite random 
variable V(b) such that sup t>0 f^3t < V(b), thus 

(5.12) -In I M t < -lnV(b) - — a.s. 

t J t 2 

A 

Taking the lim sup^^ completes the proof. □ 

6. An example: The Ornstein-Uhlenbeck process 

Here a is the identity matrix and m{x) = —x. To begin with, observe that the 
drift is generated by the potential 

(6.1) ip(x) = e 2 
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and since £(1/V> 2 ) = 0, it follows that l/ip 2 is the right invariant probability density 
(after, of course, a normalization). Furthermore 

(6.2) V-m = -d, 
so A c = d. Now, simple computations show that 

(6.3) £ a V 2 = 0, 

(6.4) £ a V = (W- V^-^V-m) , 



(6.5) 



W • ^c) 



Because the first term on the RHS of 16.41 is non-positive, using 16.21 we conclude 
£$tp < — Next choose u = ip and check that the conditions in theorems 15 . 1 1 and 
15.21 are satisfied. Hence, 
(6.6) 



Leb{A) 



1 f 1 /* d 

< lim inf — In / ip < lim sup — In / ip < — a.s. 

t^oo t J t-*oo t J 2 

X t (A) X t (A) 

7. Conclusions 



This paper has studied the behavior of a stochastic flow generated by a posi- 
tive recurrent diffusion in M. d . Even though the main theorems are proved under 
"unattractive" conditions on the coefficients, the asymptotic focusing of $Xt(A) u ^ s 
likely guaranteed under weaker assumptions. There is evidence that only an in- 
tegrable invariant measure is needed to make Theorems 14.31 15. H 15.21 hold. The 
uniform Lipschitz assumption on the drift is purely technical: it guarantees that the 
explosion time t(x,b) of Xt verifies 

(7.1) P[e(x, b) = oo, for all x G M d ] = 1, 

which is slighty stronger than bare positive recurrence. Similarly, the reversibility 
condition, critical to prove transience of Xf, should not be necessary. Unfortunately, 
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as R. Pinsky points out in the analysis of a non-symmetric generator in R d i 



is 



quite difficult. 

The key property to exploit in future investigations seems to be the critical nature 
of the operator £" which offers a broad set of superharmonic functions. 

The challenged reader may be interested in the following problem: Let A be a 
compact set in M. d and let Xt be a reversible recurrent diffusion in M. d with smooth 
coefficients and integrable invariant density tp~ 2 . Prove that 
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